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METHOD AND DEVICE FOR RECONSTRUCTION OF A DYNAMIC 
THREE-DIMENSIONAL IMAGE OF AN OBJECT COVERED BY 
A CONTRAST MEDR7M 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of a priority under 35 USC 

119 to French Patent Application No. 0010211 filed August 2, 2000, the entire 
contents of which are incorporated by reference. 

BACKGROUND OF THE INVENTION 

[0002] The invention concerns the reconstruction of a three- 
dimensional image of an object covered by a contrast medium. 

[0003] Its application is of particular interest in the medical field, in 

which reconstruction of the internal structures of patients under examination is 
undertaken and, in particular, the reconstruction of angiographic images, that is, 
obtaining images of opacified vasculatures by injection of a contrast medium, 

[0004] The invention can, nevertheless, have applications in other 

fields, notably, in nondestructive industrial control, in which examinations of the 
same type as medical examinations are performed. 

[0005] In the medical field, two-dimensional projected images of the 

object, for example, a patient's head, are generally obtained by rotation of an X-ray 
camera turning around the object. There are essentially two types of reconstruction 
algorithms in X-ray imaging. A first type provides for a calculation by back 
projection and filtering or even a Fourier transform reconstruction in several 
dimensions. A second type concerns the iterative methods of reconstruction, also 
called algebraic. The principle of such an algebraic algorithm is disclosed, for 
example, in French Patent Applications Nos. 8903606, 8916906 or 9807371, which 
describes an application of an iterative algorithm of algebraic reconstruction of 
images on a multi-resolution volume. 
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[0006] In short, after a calibration of the apparatus used to determine, 

notably, the parameters of projection in the projection planes of the acquired images, 
of an observed volume broken down into elementary volume elements or voxels 
(those calibration parameters forming projection matrices), the algebraic image 
reconstruction algorithm is used to reconstruct the three-dimensional volume from 
those two-dimensional projected images. The basic principle of that algorithm is to 
initialize the voxels of the volume to a predetermined initial value, for example, a 
zero value, and to iterate a number of times the following operations: projection of 
voxels in the plane of each acquired image so as to obtain a virtual image, 
determination of the difference between the projected volume (virtual image) and the 
corresponding acquired image, and then back projection of that difference in volume. 
After a number of iterations, an estimated value representative of the density of 
contrast medium injected in the vessels X-rayed is obtained for each voxel, which 
makes it possible to visualize in three dimensions the cartography of those X-rayed 
vessels. 

[0007] Those three-dimensional images are of valuable assistance to 

the neurologist and the surgeon, whether for a diagnosis, planning of therapeutic 
procedures or evaluation of the shape and size of objects. 

[0008] On the other hand, such reconstructed volume images have 

one major dis-advantage. Actually, they do not make it possible to visualize 
propagation of the contrast medium injected in the arteries, since the acquired two- 
dimensional images, from which the three-dimensional images has been 
reconstructed, are acquired for a quasi-stationary condition of the contrast medium. 
In fact, for the reconstruction of a three-dimensional image, which will be described 
here under the term "static" by reason of the quasi-stationary character of the 
contrast medium, it is sought rather to obtain a set of images corresponding to a same 
degree of propagation of the contrast medium. 

[0009] Consequently, with such a reconstruction, the patient's 

vascular system cannot be analyzed with just three-dimensional information and 
blood flow information (propagation of contrast medium). 
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BRIEF DESCRIPTION OF THE INVENTION 

[0010] An embodiment of the invention is directed to reconstruction 

of a three-dimensional image that will be described here as "dynamic," in contrast to 
the static three-dimensional image mentioned above, making it possible to visualize 
in three dimensions the propagation of the contrast medium in the object. 

[001 1] An embodiment of the invention is directed to a method of 

reconstruction of a dynamic three-dimensional image of an object covered by a 
contrast medium, the method comprising the steps of: 

[0012] reconstruction (or reconstitution) of a static three-dimensional 
image of the object from a first set of digitized two-dimensional projected images of 
the object respectively obtained for different imaging positions around the object; 

[0013] acquisition of at least a second set of n static two-dimensional 
projected images respectively obtained for a same imaging position and at n 
successive propagation times of the contrast medium, and 

[0014] reconstruction of the dynamic three-dimensional image of the 
object from each static two-dimensional image of the second set and the 
reconstructed (or reconstituted) static three-dimensional image. 

[0015] The invention also proposes a device for reconstruction of a 

dynamic three-dimensional image of an object covered by a contrast medium, 
comprising: 

[0016] first means for reconstructing a static three-dimensional image 
of the object from a first set of digitized two-dimensional projected images of the 
object, respectively obtained for different imaging positions around the object; 

[0017] means for acquiring at least a second set of n static two- 
dimensional projected images from a same first position of the camera and 
corresponding to n successive propagation times of the contrast medium, and 
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[0018] second means for reconstructing the dynamic three-dimensional 
image of the object from each static two-dimensional image of the second set and 
from the reconstructed static three-dimensional image. 

[0019] The invention also concerns a device for reconstruction of a 

dynamic three-dimensional image of an object covered by a contrast medium, 
adapted for use of the method, as defined below. 

[0020] The different means making possible the use of the method 

are software means executed by a microprocessor. The software means can be 
provided on a support, such as a read only memory and/or a diskette. The invention 
therefore also proposes a computer program comprising program code means 
employing the method, as defined below, when the program is executed in a 
processor. 

[0021] The invention further proposes a support, such as a read only 

memory or a diskette, capable of being read by a processor, and containing program 
code means that can apply the method, when the program is executed- in the 
processor. 

BRIEF DESCRIPTION OF THE DRAWINGS: 

[0022] Other advantages and characteristics of the invention will 

appear on examination of the detailed description of nonlimitative methods of use 
and embodiments and of die attached drawings, in which: 

[0023] Figure 1 schematically illustrates a set of two-dimensional 
projected images around an object; 

[0024] Figure 2 illustrates in greater detail the acquisition of the two- 
dimensional projected images; 

[0025] Figure 3 schematically illustrates a set of projected static images 
acquired from a same vantage point; 
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[0026] Figures 4 and 5 represent a flow chart of an embodiment of the 
method; and 

[0027] Figure 6 schematically illustrates a variant embodiment of the 
invention. 

DETAILED DESCRIPTION OF THE INVENTION 

[0028] According to one embodiment of the method, the latter 
comprises a calibration of the apparatus in which is elaborated a virtual volume 
surrounding the object and broken down into voxels. The reconstructed static three- 
dimensional image is then made up of estimates, called static estimates, respectively 
associated with the voxels of the virtual volume, each static estimate being 
representative for the corresponding voxel of the density of contrast medium injected 
in the object The dynamic three-dimensional image is, in fact, composed of n 
elementary three-dimensional images corresponding to n propagation times of the 
contrast medium. A current elementary three-dimensional image corresponding to a 
current propagation time is reconstructed (1) from the static estimates of density of 
contrast medium injected in the object, (2) for each voxel, from the intensity of the 
pixel of the static two-dimensional image corresponding to the current propagation 
time, on which the voxel is projected, and (3) for each of those pixels, from the 
integral of the estimates of density of the voxels of the virtual volume situated along 
the line of sight associated with that pixel. 

[0029] In some cases, several vessels can be projected on the same 

pixel. This generally corresponds to 10% of the voxels. It can result in ambiguity in 
estimation of the reconstructed dynamic three-dimensional image. Such ambiguity 
can be resolved by using another series of static two-dimensional images acquired 
from a vantage point different from that of the n two-dimensional images of the 
second set. 

[0030] In other words, according to another embodiment, provision is 
made for a phase of acquisition of a third set of two-dimensional projected images 
respectively obtained for a same second position of the camera, separate from the 
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first position (the one corres-ponding to acquisition of the two-dimensional projected 
images of the second set), and at the same successive propagation times of the 
contrast medium. The dynamic three-dimensional image of the object is then 
reconstructed from the reconstituted static three-dimensional image, from each 
image of the second set and from each image of the third set. 

[0031] Although the invention is not limited thereto, the application of 
the method is described for the reconstruction of a dynamic three-dimensional 
angiographic image of a patient and, in particular, of the patient's head. 

[0032] Referring, in particular, to Figures 1 and 2, the imaging system 
usable for applying the invention makes it possible to obtain, notably, a set of two- 
dimensional acquired images IAl-IAm obtained in this case by the rotation of an X- 
ray source 2 around a patient's head 1 . 

[0033] As is standard in angiography, each acquired image IAi is a 
subtracted image which is, for example, obtained by a standard technique of 
logarithmic subtraction of two X-rays taken at the same angle of incidence before 
and after an injection of a contrast medium in the vasculature, the three-dimensional 
image of which it is desired to reconstruct, first static and then dynamic as described. 

[0034] More precisely, two sequences are acquired by rotation around 
the patient's head before and after injection of the contrast medium. In order to do 
so, the radiography apparatus, for example, makes two rotations of 193.5° in 
approximately 10 seconds, making possible the acquisition of 44 two-dimensional 
images taken every 4.5°. That acquisition of images IAi is intended for the three- 
dimensional reconstruction of an image of the vasculature, called static, that is, in 
which propagation of the contrast medium is not concerned. Under those conditions, 
that acquisition is made preferably when the contrast medium is quasi-stationary in 
the vessels studied. 

[0035] Each acquired image IAi is obtained from a two-dimension 
radiation detector, for example, of the luminescence amplifier type used in radiology, 
placed opposite the X-ray tube in a plane called "projection plane PPi." The 
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different projection planes are obtained by the different angular position of the 
detector in rotation around the patient's head. 

[0036] The normal line XOi to the projection plane PPi defines the 
optical axis of the acquired image IAi. The detector is connected to processing 
means 3 containing, notably, sampling means connected to a microprocessor 
incorporating as software in its associated program memory the image reconstruction 
algorithm used in the invention and, in general, all of the functional means making 
possible the use of the method, such as the different image reconstruction means. 

[0037] In the case of an X-ray imaging system, consisting of an X-ray 
source with two-dimensional detector, the geometric operation entering into 
production of the acquired image is a conical projection of an analyzed object, 
deployed in three-dimensional space, on a two-dimensional space which is that of 
the projection plane corresponding to the detection plane. The geometric parameters 
completely describing the different conical projections must be known. It is often 
impossible and too imprecise to access these parameters directly, that is, for 
example, by directly measuring on the system of acquisition the distance between the 
X-ray source and the detector. 

[0038] The operation resulting in indirect knowledge of the geometric 
parameters entering into the production of an image is called "calibration" of an 
imaging system. 

[0039] The principle, standard and known, is based on the use of a 
known geometric phantom in three-dimensional space, the two-dimensional 
projection of which is acquired. More precisely, calibration entails the following 
stages: (1) a known object is used, the calibration phantom, presenting a number of 
characteristic points whose position in space is known by coordinates measured in 
relation to a mark peculiar to that object; (2) the image of that phantom is acquired 
under the geometric conditions of a vantage point (or angle of incidence) that it is 
sought to calibrate; (3) the projections of the characteristic points are recognized in 
the image. For that purpose, each characteristic point of the object is joined to its 
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trace in the acquired image projected; (4) the system of equation describing the 
projection is inverted; and (5) and the set of parameters of the projection is, finally, 
obtained for the given vantage point. 

[0040] A form of geometric calibration phantom often used is that of 

a cube, with eight corners at which metal balls opaque to X-rays are placed. 
Calibration being an operation known to the expert, it will not be described more in 
detail. 

[0041] A known method of automatic geometric calibration of an 

imaging system by X-rays can, nevertheless, be mentioned, such as that described in 
French Patent Application No. 9300804. In brief, a phantom is used for such 
automatic calibration, in which the balls are distributed from place to place in a 
sequence such as heights of balls, measured along the axis of rotation of the imaging 
system, and especially an axis of the phantom, or monotones increasing (or 
decreasing) with a serial number of balls in the sequence. 

[0042] The calibration of the imaging system makes it possible, 

notably, to determine the estimated mean axis Ax of rotation of the camera around 
the object, as well as the position of the source 2 and the geometric characteristics of 
the optical axes of the different acquired images. Calibration also makes it possible 
to define a virtual volume W (intersection of the different projection cones) 
surrounding the object 1 and broken down into elementary volume elements Vi or 
"voxels". That volume W, and, therefore, each voxel Vi, is spatially marked in a 
reference system hereinafter called calibration reference system, one of the axes of 
which, in this case axis Z, is merged with the estimated axis of rotation Ax. It is to 
be noted here that the projection planes PPi in which the acquired images IAi are 
projected are not generally parallel to axis Z. 

[0043] Calibration also enables a projection matrix Pi to be defined 

for each acquired image IAi, making it possible to determine, for each voxel Vi, the 
coordinates of its projection (pixel) in the corresponding acquired image IAi. 
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[0044] The phase of reconstruction of the three-dimensional image 
called "static" is now going to be described, with reference in particular to Figure 4, 
from the first set of projected static two-dimensional images obtained IAi. 

[0045] An algebraic image reconstruction algorithm is going to be 

applied directly on the images IAi (obtained after calibration 30 and acquisition 31), 
so as to deliver the reconstituted static three-dimensional image IF, typically after 
three iterations. 

[0046] It is recalled that the basic principle of the image 

reconstruction algorithm used here consists of initializing the voxels of the volume at 
a predetermined initial value, for example, zero value, and then performing a number 
of iterations, typically three. Each iteration involves, for each acquired image, the 
projection of the voxel in the plane of each acquired image, so as to obtain a virtual 
image, determination of the difference between the projected volume (virtual image) 
and the corresponding acquired image, and then the back projection of that 
difference in the volume. 

[0047] Generally, after three iterations, one obtains for each voxel of 

coordinates (ij,k), an estimate ci,j,k, designated here by the term static estimate, 
representative of the density of contrast medium injected in the vessels X-rayed, 
which makes it possible to visualize in three dimensions the cartography of those X- 
rayed vessels (image IF). 

[0048] Of course, within the scope of this invention, it is not 

necessary to visualize the static three-dimensional image EF thus obtained on the 
display screen of the imaging system, since, in fact, it consists here only of an 
intermediate stage for obtaining, as is now going to be seen more in detail, the 
reconstructed dynamic three-dimensional image. 

[0049] The expert may refer, if so desired, for more details 

concerning the use of such an algebraic algorithm of image reconstruction, to the 
aforementioned French Patent Application No. 9807371. 
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[0050] In addition to acquisition of the first set of images IAi, the 

acquisition of a second set of static two-directional projected images IBl-IBn (Figure 
3) is undertaken. That series of images IBi is acquired for a same position of the 
camera, corresponding to a same vantage point PDV1. 

[0051] The n images IBt thus acquired correspond to n successive 

propagation times of the contrast medium. 

[0052] As illustrated in Figure 5, from those images IBi acquired in 

the course of an acquisition stage 50, and from the reconstituted static three- 
dimensional image IF, the three-dimensional dynamic image of the vasculature is 
going to be reconstituted in a stage 5 1 . 

[0053] It is to be noted here that the order in which the acquisitions 

31 (Figure 4) and 50 (Figure 5) are made after calibration does not matter for use of 
the method according to the invention. 

[0054] The reconstituted three-dimensional dynamic image IFD is 

composed, in fact, of n elementary three-dimensional dynamic images, on each of 
which it is going to be possible to visualize the state of propagation at the given 
propagation time t of the contrast medium in the vasculature. By successively 
visualizing the n elementary three-dimensional images, it is going to be possible to 
visualize dynamically the propagation of the contrast medium in the vasculature. 

[0055] Similar to what has been described above for the static 

reconstituted three-dimensional image, visualization of the dynamic reconstituted 
three-dimensional image is obtained by an estimate, for each value of coordinates (i, 
j, k> of the density of contrast medium at propagation time t That estimate, valid at 
propagation time t, is designated by ci, j,k (t). 

[0056] And that estimate is defined by formula ( 1 ) below: 

E (i) 
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[0057] in which: 

[0058] ci,j,k designates the static estimate of contrast medium of that 

same voxel in the static reconstituted three-dimensional image IF, 

[0059] p(t) designates the intensity of the pixel of the image IBt 

acquired at the propagation time t and on which the voxel of coordinates i j,k is 
projected, and 

[0060] E designates for that pixel p the integral of the static estimates 

of density of the voxels of the virtual volume W situated along the line of sight 
associated with that pixel p. 

[006 1 ] E is defined by formula (2) below: 

[0062] E = Scij,k Xij,k i„j,k (2) 

[0063] in which X i,j,k equals 0 if the voxel of coordinates i,j,k is not 
situated on the line of sight and equalsl if it is situated on the line of sight. 

[0064] It is to be noted here that determination of the pixel p is 

possible by reason of knowledge of the projection matrix associated with the series 
of static images taken along a predetermined line of sight. Likewise, the values A. in 
formula (2) are also known from the values of that projection matrix. 

[0065] It is to be noted here that if at time t the voxel of coordinates 

i,j,k is not yet reached by the contrast medium, then there must be a low value for 
ci, j,k(t). If value p(t) is low, that is, if the density of the pixel projected is low, 
formula (1) is therefore going to well result in a low voxel density at time t. 

[0066] On the other hand, if value p(t) is high, that means that the 

line of sight concerned crosses other luminous voxels which come to parasite the 
density due to the voxel of coordinates ij,k concerned. Formula (1) then renders the 
voxel more accessible than it really is. 
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[0067] If the voxel of coordinates i,j,k is reached by the contrast 

medium, formula (1) can create artifacts, as in the preceding case, if the voxel 
considered is superposed with a voxel not yet reached by the contrast medium. 

[0068] However, the above formula (1) makes possible a very good 

estimate of the dynamic reconstituted image, considering that the conditions of 
creation of artifacts mentioned above concern only approximately 10% of the voxels 
of the volume. 

[0069] Nevertheless, in order to resolve that ambiguity and 

consequently improve dynamic reconstruction of the image, one can, as mentioned 
in Figure 6, undertake the acquisition of another series of n static two-dimensional 
images projected here from a vantage point different from the vantage point for the 
first series of images IBi. Images IBi ICi are respectively acquired at the same 
times . 

[0070] In that case, as indicated in Figure 5 in dotted lines, 

reconstitution of the three-dimensional dynamic image in stage 5 1 is also carried out 
from the images ICi acquired in stage 52. 

[0071] A first value ci,j,k(t) will then be calculated for voxel i,j,k at 

time t, using images IBi, and a second value ci,j,k(t) will be calculated by this time 
using formula (1) in together with images ICi. 

[0072] And the final value ci,j,k(t) associated with the voxel of 

coordinates i, j,k will then be equal to the lower of the two values mentioned above. 

[0073] Various modifications in structure and/or steps and/or function 
may be made by one skilled in the art without departing from the scope and extent of 
the invention as recited in the claims. 
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